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ABSTRACT 


A rational  phenomenological  theory  of  fatigue  lifetime  prediction  under 
arbitrary  variation  of  cycle  amplitude  is  presented.  The  theory  is  based  on  the 
concept  of  damage  curve  families  and  on  an  equivalent  loading  postulate  which  de- 
fines specimens  that  have  suffered  identical  damage  under  different  loading  pro- 
grams. Lifetime  analysis  has  been  performed  for  various  cases  of  piecewise  con- 
stant and  continuous  variation  of  cycle  amplitude.  For  continuous  variation,  the 
method  requires  numerical  integration  of  nonlinear  first  order  differential  equations 
which  have  been  established  in  this  work. 
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1.  Introduction 

A basic  problem  in  design  for  fatigue  safety  is  the  prediction  of  lifetime  of 
a structural  part  when  the  amplitude  of  cyclic  load  varies  in  prescribed  fashion 
with  the  number  of  cycles. 

The  classical  test  for  fatigue  failure  is  sinusoidal  constant  amplitude 
cycling-to-f ailure  of  a specimen  or  structural  part  resulting  in  the  number  of 
cycles-to-failure  or  lifetime  N(o),  where  0 is  the  stress  amplitude.  The  plot 
of  a against  N is  known  as  the  S-N  curve.  Suppose  now  that  the  specimen 
is  subjected  to  a loading  program  cr(n)  where  n is  number  of  cycles.  The 
problem  is  to  determine  N - the  number  of  cycles  to  failure  under  this  loading 
program. 

In  view  of  the  extreme  complexity  of  internal  fatigue  failure  mechanisms, 

there  is  little  hope  of  resolving  the  problem  on  the  basis  of  micro-structural 
\ * 
considerations.  An  alternative  is  to  consider  the  problem  in  phenomenological 

fashion.  It  should,  however,  be  borne  in  mind  that  fatigue  failure  test  data  are 

subject  to  considerable  scatter.  In  order  not  to  further  complicate  a vei  ’ diffic  ilt 

problem,  it  is  reasonable  to  disregard  initially  the  statistical  scatter  aspects  o 

•i 

the  problem.  Thus  all  results  are  to  be  interpreted  in  some  "central  measure"  % 

sense  The  problem  considered  in  this  sense  has  become  known  as  cumulative  damage  i 

-1 

theory . It  has  been  the  subject  of  numerous  investigations,  particularly  in  the 
last  thirty  years.  Since  there  are  many  excellent  reviews  of  previous  work,  e.g. 

1 1,  2,  3,  4],  present  discussion  will  be  limited  to  some  selected  aspects. 

A basic  concept  in  the  approaches  given  is  the  damage  function  which  defines  j 

n some  sense  the  damage  produced  in  a specimen  when  subjected  to  n cycles  at 
stress  amplitude  O . This  function  is  often  written  as  F(n/N)  where  N is  the 


lfetime  for  constant  amplitude  O cvclinc.  Past  work  seems  to  ha\re  K«>«n  limited 
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to  the  case  where  the  loading  is  piecewise  constant,  here  referred  to  as  multi- 
stage loading.  Particular  attention  has  been  given  to  the  simplest  case  of  two- 
stage  loading.  For  multi-stage  loadings  the  damage  function  is  written 


F( 


(1.1) 


It  is  stipulated  that 


0 < F _<  1 


(1.2) 


tailure  being  defined  by  F=l. 

A simplistic  and  well  known  damage  function  has  been  postulated  by  Palmgren 
[5]  and  by  Miner  [6]  and  has  become  known  as  Miner's  rule.  According  to  them, 
failure  in  a multi-stage  loading  is  defined  by 

l Jr = 1 (1-3) 

i I 


vhere 


n^  - number  of  cycles  at  level  0^ 

N.  - lifetime  at  constant  amplitude 

According  to  this  rule  damage  produced  by  0^  for  n^  cycles  is  defined  as 
a^/N  and  the  individual  damages  are  additive  and  independent  of  sequence. 

This  simplistic  approach  does  not  in  general  comply  with  reality.  It  is  known 
( xperlmentally  and  it  is  physically  plausible  that  the  order  of  loading  can  signif- 
icantly affect  the  lifetime  of  a specimen.  Thus,  in  two-stage  constant  amplitude 
loadings,  the  left  side  of  (1.3),  which  is  sometimes  termed  Miner's  coefficient. 
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. s usually  different  from  unity  and  is  dependent  on  loading  sequence. 


It  should  be  noted  that  the  only  material  characteristic  entering  into  Miner's 
■ule  is  the  S-N  curve  for  constant  amplitude  loading.  It  seems  unreasonable  to 
assume  that  such  simple  information  would  be  sufficient  for  lifetime  prediction 
under  arbitrary  loading  programs. 

Marco  and  Starkey  [7]  assumed  that  "damage"  produced  by  n cycles  at  level  a 
ct 

s given  by  (n/N)  where  the  exponent  is  a function  of  CT.  They  arrived  at  a 
failure  condition  which  for  two-stage  loading  assumes  the  form 


n a /ot-  no 

( — ) 1 + — = 1 
v N N 

1 2 


(1.4) 


■ nd  can  be  generalized  to  multi-stage  loading.  Apparently,  the  dependence  of  the 
exponent  a on  0 must  be  uncovered  by  experimental  means  and  this  does  not 
appear  to  be  an  easy  task.  It  will  be  seen  that  (1.4)  is  a specific  special  case 
of  the  theory  to  be  developed  in  this  work  which  does  specify  the  dependence  of 
u on  0 in  unequivocal  fashion. 

Shanley  ['ll  has  constructed  a theory  in  which  it  is  postulated  that  damage 

Ql 

an  be  described  in  terms  of  Co  n.  Corten  and  Dolan  [9]  described  damage  in  the 
form  Cn  . A detailed  account  of  these  and  other  approaches  may  be  found  in  [1]. 

The  present  work  provides  a systematic  new  approach  to  the  problem  which  is 
based  on  the  concept  of  damage  curve  families,  to  be  defined  further  on.  The 
cescription  of  "damage"  by  some  function  is  avoided.  A better  concept  is  residual 
lifetime,  this  being  a quantity  which,  unlike  damage,  can  be  defined  and  measured. 
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The  damage  curves  are  defined  in  terms  of  this  quantity. 
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. The  Concept  of  Damage  Curves 

We  consider  a specimen  which  is  subjected  to  cyclic  loading  of  constant 
maximum  amplitude  0^ , failure  taking  place  after  N^  cycles.  The  pair  of  numbers 
II  ; 0^  define  a point  on  the  S-N  curve  of  the  specimen,  Fig.  1. 

If  the  specimen  is  subjected  to  n^  cycles,  n^  < N^  , at  , it  suffers  a 
i ertain  amount  of  unknown  internal  damage.  The  pair  of  numbers,  n^  ; 0^  , defines 
a point  in  the  s-n  plane.  The  region  bounded  by  the  n,0  axes  and  the  S-N  curve 
nay  be  considered  as  a damage  region.  Each  of  the  points  within  the  damage  region 
specifies  in  some  sense  an  amount  of  damage  for  n cycles  at  amplitude  0 . 

Next  we  consider  the  situation  where  after  n^  cycles  at  > cycling  is 
ontinued  at  some  stress  amplitude  0.  The  remaining  number  of  cycles  to  failure 
s called  the  residual  lifetime  and  is  denoted  N-n  where  N is  the  constant 
amplitude  lifetime  at  0.  The  pairs  of  numbers  n;0  define  a curve  in  the  S-N 
plane,  Fig.  1,  which  is  called  the  damage  curve . Evidently,  the  damage  curve  undei 
consideration  must  pass  through  the  point  ; 0^. 

It  should  be  noted  that  in  this  description  the  damage  suffered  by  a specimen 
s essentially  defined  by  its  residual  lifetime  which  is  a measurable  quantity, 
he  S-N  curve  is  the  ultimate  damage  curve  when  remaining  lifetime  is  zero. 

The  basic  assumption  is  made  that  a damage  curve  is  uniquely  defined  by  its 
nitia]  point.  This  implies  that  if  any  other  point  on  the  damage  curve  n^ : , 

■ay,  is  interpreted  as  the  first  stage  of  a two-stage  cycling  experiment,  the 
emaming  lifetimes  for  various  subsequent  stages  will  be  determined  by  the  initia 
carnage  curve  which  passes  through  n^;  . This  assumption  can  be  accepted  with  or 

without  distinction  between  low-high  and  high-low  loadings.  Further  discussion  of 


i h is 


latter  aspect  is  given  later. 


* 
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Uniqueness  of  the  damage  curve  passing  through  a point  and  subsequent  analytical 
procedures  for  lifetime  computations  are  based  on  a fundamental  postulate  which  will 
now  be  explained.  For  this  purpose  we  shall  introduce  the  concept  of  equivalent  cyclic 
loadings . Cyclic  loading  programs  are  termed  equivalent  for  stress  amplitude  0 if 
for  subsequent  constant  amplitude  loading,  at  same  stress  level  0,  residual  lifetimes 
are  the  same.  The  equivalent  loading  postulate  is  now  stated  as  follows:  Cyclic  load- 

ings which  are  equivalent  for  one  stress  level  are  equivalent  for  all  stress  levels. 

To  show  uniqueness  of  damage  curves  on  the  basis  of  this  assumption,  consider 
the  damage  curve.  Fig.  1,  as  having  been  obtained  by  a series  of  two-stage  loadings 
in  all  of  which  the  first  stage  is  n^,  0^,  as  previously  explained.  In  the  two 

stage  program  for  n^  cycles,  02  for  an  infinitesimal  number  of  cycles  the  resi- 

- r r 

dual  lifetime  is  n2  . In  the  loading  program  O2  for  Nj  - ^ cycles 'the  residual 

F 

lifetime  is  evidently  also  ^ . Therefore  these  two  loadings  are  equivalent.  It 
follows  by  the  equivalent  loading  postulate  that  if  after  completion  of  these  two 
loading  programs  the  specimens  are  cycled  to  failure  at  some  other  stress  level,  o, 

F 

the  residual  lifetime  will  be  the  same  n , Fig.  1.  Therefore  the  loadings 

F 

for  n^  cycles,  a for  infinitesimal  number  of  cycles;  f°r  ^ cycles,  a for 

F 

infinitesimal  number  of  cycles  have  the  same  residual  lifetime  n for  subsequent 
constant  a cycling  to  failure.  This  proves  uniqueness  of  the  damage  curve. 

A damage  curve  as  defined  describes  the  results  of  a large  number  of  fatigue 

tailure  experiments  in  two  level  loadings.  It  is  desirable  that  a maximum  number 

of  parameters  which  influence  fatigue  failure  be  kept  constant  in  this  series  of 

experiments.  The  parameters  to  be  considered  are  : frequency  of  loading,  mean 

stress  0 - the  average  of  maximum  and  minimum  amplitudes  in  a cycle  and  R - 

m 

the  ration  between  minimum  and  maximum  amplitudes  in  a cycle.  It  may  be  assumed 
that  the  frequency  is  constant  in  all  cyclic  ioadings  and  that  either  or  R is 
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constant.  In  the  special  case  of  a ■*  0,  R “ -1  both  0 and  R can  be  kept 

in  m 1 

constant.  It  is  In  general,  to  be  expected  that  a change  in  any  of  these  parameters 
will  change  the  damage  curves. 

Because  of  the  inherent  scatter  of  fatigue  testing,  all  damage  curves  must 
be  interpreted  in  some  best  fit  sense,  as  is  the  case  for  S-N  curves. 

The  equations  of  the  damage  curves  are  written  in  the  general  form 


o = a ( n , y ) = asf(n,y) 


(2.1) 


where  Y is  a parameter,  0g  is  the  static  failure  stress  or  some  other  con- 
venient stress  parameter  and  f is  a non-dimensional  function.  We  now  proceed 
to  discuss  some  general  properties  of  the  damage  curves. 

(.a)  All  damage  curves  pass  through  the  point  0;  0^ . __ 

Evidently,  since  static  failure  occurs  at  one  quarter  cycle  which  is  considered 
as  n~0,  it  follows  from  (2.1)  that 


KO, Y)  = 1 


(2.2) 


However,  curve  fitting  difficulties  mav  require  use  of  a a value  which  is  not 

s 

equal  to  the  static  test  result. 

(b)  Damage  curves  do  not  Intersect. 

To  prove  this  assertion  we  consider  the  case  of  two  low-high  load  cycling  pro- 
> rams  with  equal  amplitude  pairs  and  a The  number  of  cycles  at  0^  is  n^ 

in  the  first  program  and  ' > n^  in  the  second  program.  Obviously  the  remaining 
I He  time  - n^  at  0^’  in  the  first  program  must  be  larger  than  the  remaining 
lifetime  Nj  - n^'  in  the  second  program.  Consequently, 


n?'  > n2  (2.3) 

Pig  2a  shows  a hypothetical  intersection  between  two  damage  curves.  an<^ 

IS  B„  represent  the  two  loading  programs  described  above.  It  is  seen  that  (2.3) 

i 2 

i , iolated.  Therefore,  intersection  is  not  possible  and  neighboring  damage  curves 
re  as  in  Fig.  2b. 

. ) Damage  curves  do  not  intersect  the  n or  a axis  (except  at  a = a ) . 
s 

Suppose  that  a damage  curve  intersects  the  n axis.  This  wou . imply  that 

n a two-stage  loading  program  with  amplitudes  0,0  the  lifetime  is  smaller  than 

’J(.C).  But  the  lifetime  must  be  N(o)  since  the  loading  is  teally  one  stage. 

herefore  the  damage  curve  must  pass  through  the  point  0,0  which  is  not  possibl  i 

m view  of  (a) . 

Suppose  that  a damage  curve  passing  through  n^,  0^  intersects  the  C axis 
at  a . This  would  imply  that  if  the  specimen  is  exposed  to  n^  cycles  at 
ne  remaining  lifetime  for  cycling  at  o^  is  N^,  which  can  only  be  correct  if 
: | =0.  Therefore,  no  damage  curve  intersects  the  0 axis  below  o . 

in  view  of  the  properties  (a)  - (c)  of  the  damage  curves,  it  is  seen  that  th >y 

rm  a family  of  non- intersecting  curves  within  the  region  of  the  s-n  plane 

hich  is  bounded  by  the  n,a  axes  and  by  the  S-N  curve.  The  damage  curves  ail 

kiss  through  the  point  0;O  and  approach  the  n,0  axes  asymptotically.  The  S-N 

s 

urve  is  the  limiting  damage  curve.  It  corresponds  to  a two  stage  loading  in  whi  h 
atlure  occurs  for  the  first  stress  and  residual  lifetime  for  0 is  zt  r . 

The  simplest  procedure  to  assess  the  form  of  the  damage  curves  is  to  assume 
hat  their  mathematical  description  is  similar  to  that  of  the  S-N  curve.  Let 


'he  equation  of  the  S-N  curve  no 


represented  in  the  form 
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o = o f(N,F) 
s 


(2.4) 


here  N denotes  value  of  n for  a point  on  the  S-N  curve.  Then  the  equation 
of  the  damage  curves  is  assumed  to  be  given  by  the  similar  form 


o = Ogf(n,y)  (2.5) 

The  parameter  y for  a damage  curve  passing  through  the  point  n^;C^  is  determini  d 
by  inserting  those  values  into  (2.5)  and  solving  for  y. 

The  two  most  common  analytical  representations  of  S-N  curves  are: 

ti)  O is  a linear  function  of  log^N.  This  will  be  called  the  semi-log  form. 

(2)  log^o  is  a linear  function  of  log^N.  This  will  be  called  the  log-log 
form.  Thus 


0 = a (1  + T log  N) 


(2.6) 


or 


0 = 


0 N1 
s 


In  both  of  these  representations  N=1  is  considered  physically  equivalent  to 
N-0 

It  follows  that  the  corresponding  forms  of  damage  curves  passing  through  a 
specified  point  n^ ; 0^  are: 


o = a (1  + y n) 


(a) 


a Jo  -1 

1 s 
log  n 


(b) 


(2.8) 
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o = o n 
s 


log(a/o  ) = y log  n 


(a) 


log(a1/os) 

log  n1 


(b) 


(2.9) 


The  form  (2.8)  is  generally  more  suitable  for  fiber  composites  while  the  form 
(2.9)  is  more  suitable  for  metals. 

The  analytical  expressions  (2. 6-2. 9)  are  merely  convenient  approximations  and 
should  not  be  taken  too  literally.  It  is  seen  that  all  curves  intersect  the  N 
axis  contrary  to  theorem  (c).  The  conclusion  is  that  for  very  low  stress  ampli- 
tudes the  analytical  expressions  are  not  good  approximations. 

Other  fitting  difficulties  may  occur  at  small  number  of  cycles.  It  fre- 
quently happens  that  the  determined  from  S-N  curve  fitting  is  not  the 

0 found  bv  static  testing.  In  that  case  the  former  should  be  chosen  and  be 
regarded  as  a curve  fitting  parameter. 

It  is  permissible  to  fit  a convenient  portion  of  the  S-N  curve,  for  the 
range  i o I 0^,  say,  by  the  analytical  expressions.  In  that  case,  the  valid- 
ity of  the  damage  curve  expressions  is  limited  to  that  loading  range  and  is 

a fitting  parameter. 

The  S-N  curve  may  be  represented  by  more  than  one  function.  For  example: 
a broken  straight  line  in  log-log  presentation.  In  that  case  the  damage  curves 
would  be  similarly  represented. 

If  the  analytical  forms  do  not  prove  accurate  enough  it  is  necessary  to  deter- 
mine the  damage  curves  by  testing.  In  order  to  determine  a damage  curve  which  passes 
through  the  point  n^c^  3 RrouP  specimens  are  subjected  to  n^  cycles  at 
and  are  then  cycled  to  failure  at  different  stress  levels  0^,  0^,  ...  . Rach  of 
tiie  lifetimes  at  0^,0^,  ...  determine  a point  on  the  damage  curve. 
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l'he  damage  curves  can  be  separately  defined  for  two  stage  loadings  which 
are  of  the  low-high  type  or  of  the  high-low  type.  There  then  arises  the  question 
if  both  sets  of  curves  are  identical.  Consider  for  this  purpose  a loading  pro- 
gram in  which  n^  cycles  take  place  at  o^.  Then  the  amplitude  is  raised  to 
v for  the  duration  of  a very  small  number  of  cycles.  After  that,  the  amplitude 
is  reduced  again  to  ci  ^ , and  cycling  is  continued  until  failure.  A diagram  of 
such  a loading  program  in  the  s-n  plane  is  shown  in  Fig.  3.  If  the  number  of 
cycles  to  failure  is  insignificantly  affected  by  the  rise  to  o^,  and  is  thus 
equal  to  N|,  the  lifetime  at  then  the  segment  AjA7  is  insignificant  and 

tile  upward  damage  curve  A ^ H ^ coincides  with  the  downward  damage  A7B^.  If, 
th<wever,  this  is  not  the  case,  the  upward  and  downward  damage  curves  are  ditferent 
and  the  damage  curves  exhibit  a hysteresis  phenomenon.  On  the  basis  of  experience 
in  tutigue  testing,  it  is  reasonable  to  assume  that  such  an  effect  will  not  be 
significant  at  low  to  moderate  stress  rises  but  may  be  significant  at  high 
stress  lises  with  consequent  conclusions  for  the  identity  of  upward  and  downward 
damage  curves. 

It  is  to  be  emphasized  that  it  does  not  seem  possible  to  include  such 
hysteresis  effects  in  the  simple  analytic  representations  (2.8  - 2.9)  since  they 
.re  bused  on  the  form  of  the  S-N  curve  where  hysteresis  cannot  enter.  There- 
•re,  such  effects  must  be  uncovered  by  experiments.  This  appears  to  be  a 
difficult*  task  because  of  the  considerable  scatter  of  fatigue  failure  results 


11 


3 . Analysis  of  Loading  Programs 

The  simplest  loading  program  is  two-stage  cyclic  loading.  Its  analysis 
in  terms  of  the  damage  curves  is  trivial  since  the  damage  curves  have  been 
defined  in  terms  of  residual  lifetime  under  such  loadings.  Let  0^  be  applied 
for  n^  cycles.  Then  the  stress  amplitude  is  changed  to  O ^ and  cycling  is 
continued  until  failure  occurs  after  n ^ cycles.  It  is  required  to  compute  n^ , 

Let  it  be  assumed  that  the  S-N  curve  is  of  form  (2.6)  and  that  the  damage 
curves  are  of  semi-log  form  (2.8a).  Application  of  n^  cycles  at  o ^ defines  a 
damage  point,  Fig.  l.  The  damage  curve  passing  through  this  point  is  given  by 
(2.8).  It  follows  that 


log(N„-n„)  O 

0o  = 0 [1  + — — — ( — - 1)] 

2 sl  log  n,  a ' 

I s 


(3.1) 


Solving  for  we  have 


n2  = N2  " "I1""1 


s,  = — 
1 o 


S2  o 


(3.3) 


[’he  quantity 


n n 

M = — + — 
N1  N2 


(3.4) 


is  defined  as  Miner's  coefficient  (for  two  stage  loading).  It  follows  easily 


from  (2.6a)  and  (2,8a)  that  (3.2)  can  be  written  as 


which  replaces  Miner's  rule  for  the  present  case.  It  follows  from  (3.4  - 5)  that 


1-s- 


M = 1 + 


- ( 


-I  \1-S1 


-t  is  easily  shown  that 


(3.6 


M > 1 


S1  * S2 


when 


(3.7 


M < 1 


S1  * S2 


A similar  analysis  for  damage  curves  of  log-log  type,  (2. 9), yields  the 
resul  t s 


log  s2/log  sx 

1 ”7*  + 


n,  n.  108  S2/l0g  S1 

M=  1+^  - 


(3.8 


Again,  it  is  easily  shown  that  (3.8)  complies  with  (3.7). 

Note  that  (3.5)  and  (3.8)  are  of  the  form  (1.4)  which  is  based  on  an  expon- 
ential damage  rule  postulated  by  Marco  and  Starkey  [7], 

Plots  of  Miner's  coefficient  for  two-stage  loading  as  given  by  (3.6)  and 
(3.8)  are  shown  in  Fig.  4.  It  is  seen  that  there  is  considerable  difference 
between  the  present  prediction  and  Miner's  rule. 
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We  now  proceed  to  the  case  of  multi-stage  loadings.  The  approach  will  be 


I 

L 


illustrated  by  treatment  of  a three-stage  loading.  A specimen  is  subjected  to  n 


cycles  at  amplitude  O^,  then  to  n^  cycles  at  amplitude  o^.  Then  the  stress 


level  is  changed  to  and  it  is  required  to  determine  the  remaining  lifetime 


n^.  Consider  the  two  stage  loading  cr^,  cr^ . Once  the  stress  level  is  changed  to 


0 ^ we  proceed  on  the  damage  curve  passing  through  point  n^;  , Fig.  5,  until  the 


level  0^.  At  this  stage  the  remaining  lifetime  is  - n^*  After  n^  cycles 


are  applied  at  this  stress  level  the  remaining  lifetime  is  - n^  - n^.  It 


follows  that  if  the  specimen  is  subjected  to  n.^  + n9  cycles  at  0 the  remain- 
ing lifetime  is  the  same  as  that  of  a specimen  subjected  to  the  two  stage  loading 
under  consideration.  (In  much  vaguer  language  we  would  say  that  the  two  loadings 
result  in  the  same  damage).  Consequently,  the  three  stage  loading  is  equivalent 


to  the  two  stage  loading  a ^ n^  + n2  cycles,  -*•  n^  cycles.  On  the  basis  of 


the  equivalent  loading  postulate  we  can  use  the  two  stage  loading  procedure  which 
implies  that  after  n ^ cycles  at  C we  proceed  to  level  on  the  damage  curve 

which  passes  through  the  point  n^  + n^ ; . This  determines  the  remaining  life- 


time n^. 


Evidently  this  procedure  can  be  generalized  to  multi-stage  loadings  with  anv 
number  of  stages.  It  follows  that:  the  damage  curves,  defined  for  two-stage 

loadings,  are  sufficient  information  to  determine  lifetime  for  any  multi-stage 


loading. 


For  analytical  purposes  it  is  convenient  to  adopt  the  following  notation:  the 
damage  curves  which  are  needed  for  transition  from  one  loading  level  to  another  tr< 


denoted  , D9,  ...,  the  Last  one  being  the  S-N  curve.  The  abscissa  of  a point 


on  damage  curve  at  stress  level  o ^ is  denoted  n^,  • Lifetime  at  i.e. 


abscissas  of  points  on  the  S-N  curve  is  denoted  N . 
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Consider  a three-stage  loading.  Fig.  5,  with  damage  curves  of  semi-log 
type.  With  the  notation  adopted 


nl  ' nil 


n2  n22  “ n12 


(3.9) 


n3  = N3  - n23 


Since  the  damage  curves  have,  by  hypothesis,  the  same  functional  form  as  the 
S-N  curve,  it  follows  at  once  that  (3.5)  can  be  interpreted  as  a relation 


between  n^,  n2  where  the  damage  curve  2 now  replaces  the  S-N  curve  of  the 


previously  considered  two  level  loadings.  Therefore 


n,  1-S2 

( )^I  + — 
21  n22 


= 1 


(3.10) 


Similarly,  we  interpret  the  points  n22'°2  anc*  n23’°3  38  3 two  ^-eve^ 


loading.  Therefore,  from  (3.5) 
l-s3 


(^2  )l-‘2+n.  ! 

N2  "3 


(3.11) 


Using  the  functional  form  of  the  S-N  and  damage  curves  (2.6)  and  (2.8a),  it 
follows  from  (3.10  - 11)  that 

l-s„ 


1-s, 


n.  i.  n 

[ ( -1  )l“8l  + -1  ]1_S 
1 N ; N.,  J 


n,, 


a 


2 + 


N 


= 1 


si  o 


(3.12) 


Therefore  the  Miner's  coefficient  is  given  by 
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n,  n n £ n y — 

M ■ 1 + r +tT-  1 ( n >1-si  +sr  1 2 

1 2 1 n2 


Similar  analysis  for  the  log-log  forms  gives  the  results 
loR32/1oes1  108  V1-* 

1 + sj>  + ‘ 1 

n,  n n log  s„/log  s n log  s./log  s„ 

M=l+  — + --[(  — ) 1 i+—  1 i 1 

N1  N2  N1  N2 


As  a final  example  for  multi-stage  loadings,  we  consider  the  case  of 
periodic  two  stage  loading  where  the  period  (called  block  in  fatigue  jargon) 
is  composed  of  for  nj  cycles  followed  by  for  cycles.  A rep- 

resentation of  this  loading  for  nondimensional  stresses  is  shown  in  Fig.  6a. 
The  loading  path  in  the  s-n  plane  is  shown  in  Fig.  6b. 

Let  the  inverse  of  a damage  curve  function  s = f(y,n)  be  written 


n = g(Y,s) 


s = o/c 


An  alternative  form  is 


y = Y(s,n) 


With  the  n,  notation  adopted  above,  we  proceed  in  the  following  manner 


along  the  Loading  path  shown  in  Fig.  6a. 
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nll 

= nl 

n12 

= g(Yj  »s2) 

Y1  = Y(n1>s1) 

n . 

— p + n 

22 

12  2 

n21 

= g ( Y 9 . s x ) 

y 2 f ^n22 * S2^ 

n31 

n21  + nl 

(3.17) 


This  step-by-step  procedure  is  continued  until  the  S-N  curve  is  reached,  thus 
determining  the  lifetime  for  the  loading. 

The  functions  g and  y for  the  two  types  of  damage  curves  considered  are 

semi-log  (3.  18) 


log-log  (3.19) 


I'he  procedure  is  very  convenient  for  numerical  computation.  It  is  also 
easily  realined  that  it  can  be  used  to  analyze  any  multi-stage  loading. 

We  now  consider  the  case  when  the  loading  amplitude  varies  in  arbitrary 
.■ntinuous  fashion  with  the  number  of  cycles,  which  is  denoted  n.  Thus 


g(s,Y)  " 

1 (s,n) 
g(&, Y) 
y(s,n) 


s— 1 

- 10  ' 

s-j 
log  n 

1/7 


log  s 
log  n 


1 


a = o(  n) 


(3.20) 
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with  inverse 

n = n(a)  (3.21) 

The  loading  (3.20)  may  be  regarded  as  a multi-stage  loading  consisting  of  an 

infinite  number  of  levels  with  increments  do.  Figure  7a  shows  the  loading  curve 

(3.20)  with  initial  value  O . The  representation  of  (3.20)  in  the  s-n  plane 

o 

is  as  follows:  At  Oq  plotted  on  the  0 axis  proceed  dn  horizontally  to  the 

right.  Trace  a damage  curve  passing  through  the  point  dnjO^  and  on  it  proceed 
until  + do,  where  do  is  determined  by  dn  and  the  slope  of  the  loading 
curve  at  0 . At  the  point  thus  reached  in  the  s-n  plane,  proceed  another  dn 
horizontally  and  repeat  the  previous  process.  Fig.  7b  illustrates  this  process 
for  a typical  point  on  the  loading  curve  (3.20)  and  the  image  of  that  point  in 
the  s-n  plane.  It  is  seen  that  by  this  procedure  the  loading  curve  C is 
mapped  into  the  curve  C'  in  the  s-n  plane.  Intersection  of  curve  C'  with 
the  S-N  curve  defines  failure  under  the  loading  program  (3.20). 

This  procedure  of  construction  of  the  image  C'  of  C will  now  be  expressed 
analytically.  Let  the  equation  of  C'  be  denoted 


G(s,n)  = 0 


s = a/a 

s 


(3.22) 


The  relation  between  da  and  dn  on  this  curve  is  illustrated  in  Fig.  7b. 
The  differential  dn  is  given  by 


dn  = dn  + dn 


(3.23) 


where  dn^  is  the 
curve  which  passes 


contribution  due  to 
through  the  point  n 


movement  Irom  0 to  O + do  on  the  damage 
^;o.  Note  that  at  this  point 
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11 D = n 


(3. 24) 


Let  the  damage  curve  equation  (2.5)  be  written  in  the  form 


nD  = g(s,y) 

Since  y is  constant  on  this  curve  it  follows  that 


d%  = It  dS 


(3.25) 


(3.26) 


in  view  of  (3.24  - 25),  y of  the  damage  curve  passing  through  the  point  nD>° 
can  be  expressed  from  (3.20)  in  the  form 


Y = Y(s,n)  (3.27) 


Now  write  (3.21)  in  the  form 


n = h(s)  (3.28) 


introducing  (3.26  - 3.28)  into  (3.23)  we  obtain 

tin  = [ + h'(s)]ds  (a) 

(3.29) 

Y = Y(s,n)  (b) 


where  (b)  is  to  be  introduced  into  (a)  af ter  the  partial  differentiation  has 
been  performed. 

Eq ii. (3-29)  is  a nonlinear  ordinary  first  order  differential  equation  of 


the  type 
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dn  ~ (t)(s,n)ds 

Ils  integration  starting  out  from  an  initial  point  n 50  gives  the  equation 
of  the  curve  C',  (3.22). 

The  intersection  of  (3.22)  with  the  S-N  curve  defines  the  stress  level  0 ^ 
at  which  failure  occurs  under  the  loading  program  (3.20).  The  cycle  ordinate  n 
tn  the  S-N  plane  associated  with  the  intersection  point  does  not  define  the 
lifetime.  The  .Lifetime.  N is  defined  from  (3.20)  by  (see  Fig.  7) 

N = n(o  ) (3.3C 

u u 

This  raises  a problem  when  there  is  no  one-to-one  correspondence  between  O and 
n as  defined  by  (3.20),  which  happens  when  (3.20)  is  not  a monotonic  function. 

In  this  event  it  is  necessary  to  determine  during  the  integration  of  (3.29)  the 
images  of  the  extrema  of  (3.20)  on  C' . In  this  fashion  is  then  located  on 

the  proper  monotonic.  part  of  (3.20). 

[n  the  important  case  of  load  discontinuities  in  (3.20)  it  is  necessary 
UJ  integrate  (3.29)  until  the  image  of  the  discontinuity  in  the  s-n  plane  is 
reached.  The  jump  in  o then  proceeds  on  the  damage  curve  which  passes  through 
the  end  point  of  integration  reached,  in  the.  s-n  plane.  Then  the  integration 
is  continued  unti L the  next  discontinuity,  etc. 

The  functions  g(y,s)  and  y(s,n)  which  appear  in  (3.29)  have  been  given 
for  semi-log  and  log-log  damage  curves  by  (3.18  - 19).  Therefore,  for  these 
cases,  Che  differential  equation  (3.29)  assumes  the  forms 
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dn  = [ An  10  + h' /s) ]ds  (a)  semi-log 

(3.31) 

dn  = ^ s "lor~s  + h'(s)]ds  (b)  log-log 

A representation  of  (3.32)  in  terms  of  the  variables  n,  n can  be  obtained 
by  writing  (3.20)  in  the  form 

s - H(n)  (3.32) 


and  substituting  into  (3.21). 

It  appears  that  analytical  integration  of  equs.  (3.31)  is  not  possible  but 
the  equations  are  readily  integrated  numerically  by  use  of  the  Runge  Kutta  method, 
for  example.  For  this  purpose  it  is  convenient  to  introduce  the  variable 

n = log  n (3.33) 


Then  (3.31)  assume  the  forms 


dn  = n 

ds  s-.l 


h’  (s) 
10n£n  10 


(a) 


dn  = . n + h1  (s)  , l 

ds  s log  s Qn  An  10  ' 


(b) 


(3.3/*) 


As  an  example  we.  consider  the  loading  program  shown  in  Fig.  8a.  The  initial 

amplitude  is  given  by  s = a fa  and  then  diminishes  linearly  to  a value  s, . 

o o s J 1 

Cycling  is  then  continued  at  constant  amplitude  until  failure.  The  linear  variation 
of  cycling  amplitude  is  described  by 
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s “ so  + (8i  " so}  T 


.t.  follows  from  the  definition  (3.28)  find  from  (3.35)  that 


(3.35) 


h'  (s)  = — (3.36) 

1 ’0 

With  this  value  of  h'(s)  equs.  (3.34)  must  be  integrated  starting  out  at  the 
initial  point  n = 0;  s = s^.  Note  that  n = 0 corresponds  to  n = 1,  not  n = 0. 
This  is  an  insignificant  correction  which  is  required  because  of  the  log  repres- 
entation. 

Integration  of  the  differential  equations  gives  the  image  C'  of  the  loading 

function  in  s-n  plane.  This  image  is  shown  in  Fig,  8b  for  the  semi-log  case, 

4 

equ. (3.34a)  when  sQ  = 0.7,  nQ  = 10  . It  is  seen  that  the  image  C'  in  the  s-n 
plane  of  the  linear  loading  curve  C starts  out  horizontally  and  then  approaches 
asymptotically  a straight  line  passing  through  the  point  r|  = 0,  s = 1.0,  i.e., 
a damage  curve.  This  is  a general  characteristic  of  solutions  of  (3.34)  since 
for  large  n the  second  term  on  the  right  sides  becomes  very  small  in  comparison 
to  the  first  terms  because  of  10^  in  the  denominator.  If  the  second  term  is 
neglected  the  solution  of  the  differential  equations  are  the  damage  curves  (2.8a), 
(2,9a)  respectively.  This  is  also  evident  from  (3.23)  since  neglect  of  the.  second 
term  on  the  right  side  of  (3.34)  is  equivalent  to  neglecting  dn  in  comparison 
t o d n . 

For  linear  loading  from  s^  to  s^  the  remaining  lifetime  is  given  by  the 
horizontal  distance  between  the  point  on  C'  where  snSp  and  the  S-N  curve,  in 
s,n  coordinates.  Note  that  the  present  figure  is  in  s,log  n coordinates.  Fig.  8b 


shows  remaining  lifetimes  for  various  values  of  s^. 
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In  order  to  compare  the  present  prediction  to  those  of  Miners'  rule  it  is 
necessary  to  write  Miners'  rule  for  a continuous  loading  function.  To  do  this 
the  continuous  loading  (3.20)  is  interpreted  as  a multi-stage  loading  with  successive 
load  increments  do  and  cycle  increments  dn.  Then  (1.3)  may  be  written 


a 

r u 


dn 


N(a) 


dn/da 

N(a) 


da 


1 


(3.37) 


where  N(a)  is  lifetime  at  constant  O as  given  by  the  S-N  curve.  Equ.  (3.38) 

defines  the  failure  stress  0 . Lifetime  N is  then  given  by  (3.30). 

u u 

In  terms  of  the  nondimensional  stress  s = a/o  (3.37)  assumes  the  form 

s 


s 

u 

f 


h'  (s) 
N(s) 


ds  = 1 


s 


0 


(3.38) 


For  semi-log  and  log-log  representations  (2.6),  (2.7)  we  have  respectively 


s-1 

r 


N(s)  = 10 

. . / s i/r 

N(s)  = s 


(3.39) 


while  h' (s)  for  the  linear  loading  is  given  by  (3.36).  For  the  present  loading 
(3.38)  assumes  the  form 


S1 

f h'(s) 

J N(s) 
s0 


ds  + 


N(s1) 


= 1 


(3.40) 


which  defines  the  residual  lifetime 


n 


1 


on  the  basis  of  Miner's  rule. 
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It  follows  from  (3.40), 


(3.36)  and  (3.39)  that 


ni  = N(si}  " I^TTo 


no  r . ^soi 
Vs!  N(sl} 


is  given  by 


(a)  semi-log 


(3.41) 


r nQ  N(^) 

■ TTI  7-^1  [s0  N(^)  "  4  Sl]  (b)  log~log 

The  n predicted  by  (3.41a)  are  also  shown  in  Fig.  8b.  It  is  seen  that 
there  are  significant  differences  between  remaining  lifetimes  as  predicted  by 
Miner's  rule  and  the  present  theory. 


n^  = N (s^) 


Analysis  has  also  been  performed  for  the  case  of  pure  linear  variation  of 
cycle  amplitude  (increasing  and  decreasing)with  number  of  cycles,  until  failure. 
Results  obtained  for  various  such  cases  showed  little  numerical  difference  between 
lifetime  predictions  by  present  theory  and  by  Miner's  rule.  It  is  of  great 
importance  to  arrive  at  some  general  conclusions  with  respect  to  loading  char- 
acteristics for  which  Miner's  rule  and  the  new  theory  predict  significantly 
different  lifetimes.  It  appears  that  much  further  analytical  work  is  needed 
to  obtain  such  information. 

4 . Summary  and  Conclusion 

A rational  phenomenological  theory  to  compute  fatigue  lifetime  under  arbitrary 
cyclic  loading  programs  has  been  developed.  The  theory  is  based  on  the  concept  of 
damage  curve  families  which  represent  all  possible  two  stage  cyclic  loadings  - to  - 
failure.  The  damage  curves  may  be  obtained  experimentally.  In  the  present  work 
their  analytical  form  is  postulated  on  the  basis  of  assumed  similarity  with  S-N 


curve  analytical  representations. 


The  theory  is  based  on  a basic  assumption  which  has  been  termed  the 


equivalent  loading  postulate  which  defines  in  macroscopic  terms  specimens  which 
have  suffered  the  same  damage  under  different  loading  programs.  This  postulate 
leads  to  uniqueness  of  the  damage  curves  and  underlies  methods  of  lifetime  analysis 
in  terms  of  the  damage  curves. 

Procedures  of  lifetime  prediction  have  been  given  for  piecewise  constant  cycle 
amplitude  variation  (multi-stage  loading)  as  well  as  for  continuous  variation  of 
cycle  amplitude  with  number  of  cycles.  In  the  first  case  analysis  consists  of 
a simple  step  by  step  procedure  in  terms  of  damage  curves.  In  the  second  case 
solution  of  initial  value  problems  for  first  order  nonlinear  differential  equations 
is  required. 

Damage  curve  families  considered  in  this  work  are  straight  lines  in  log  n,  0 
or  in  log  n,  log  0 planes. 

Analytical  results  have  shown  that  for  multi-stage  loading  programs  the 
predictions  of  the  present  theory  are  considerably  different  tiian  those  based 
on  Miner's  rule.  Analysis  of  a few  cases  of  continuous  amplitude  variation  programs 
has  shown  that  in  certain  cases  present  predictions  are  significantly  different 
from  those  based  on  Miner's  rule  while  in  others  this  is  not  the  case. 

Further  analytical  work  is  needed  to  search  for  loading  characteristics  vhuh 
lead  to  significant  discrepancies  between  prediction  by  Miner's  rule  and  present 
theory.  It  is  also  of  great  importance  to  establish  an  analytical  procedure  tor 
random  loadings,  i.e.,  where  the  cycle  amplitude  is  a random  variable  and  only 
statistical  loading  information  is  available.  This  is  a difficult  problem  since 
it  requires  solution  of  nonlinear  differential  equations  with  random  inputs. 

An  experimental  program  is  now  under  way  to  examine  the  validity  of  the  theory 
This  consists  of  two  stage  loadings-to-f aiiure  of  gloss  tiber/ep>xv  specimens. 
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